INVERTIBILITY OF SUBMATRICES OF PASCAL'S MATRIX AND 
BIRKHOFF INTERPOLATION 
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Abstract. The infinite (upper triangular) Pascal matrix is T = [(^)] for < It is easy 

to see that that submatrix T(0 : n, : n) is triangular with determinant 1, hence in particular, 
it is invertible. But what about other submatrices T{r,x) for selections r = [ro, . . . , r^j] and x = 
[xo , . . . , Xjj] of the rows and columns of T? The goal of this paper is provide a necessary and 
sufficient condition for invertibility based on a connection to polynomial interpolation. In particular, 
wc generalize the theory of Birkhoff interpolation and Polya systems, and then adapt it to this 
problem. The result is simple: T(r, x) is invertible iff r < x, or equivalently, iff all diagonal entries 
are nonzero. 
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1. Introduction. The infinite and order n upper triangular Pascal matrices are 



1111 
12 3 
13 
1 



and T(0 :n,0:n) 



1 1 
1 



1 







for i,j = 0, 1, . . ., with (-J) := if i > j. Note that r(0 : n,0 : n) can be viewed as a 
submatrix, or truncation, of T determined by choosing the first n+1 rows and columns. 
General truncations are determined by selecting arbitrary (finite) combinations of the 
rows and columns of T. These can be represented as 



T(r,x) 



for some selections r 
respectively. 



r") 



KroJ 
(xi 



O CD 



[ro, . . . , rd\ and x = [xq, . . . , Xd] of the rows and columns of T, 
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The main issue taken up in this work is the invertibility of truncated Pascal 
matrices. It is trivial to see that Tn is invertible, with determinant 1. But what 
about the others? For example, consider 



r([0,l,2],[l,2,5]) 



11 1' 
1 2 5 
1 10 



and r([l,3,4],[l,2,5]) 



12 5' 
10 
5 



The second matrix, a triangular matrix with a zero on the main diagonal, clearly has 
a zero determinant, and so is not invertible. However, the first matrix is invertible, 
but is not obvious to see. More precisely, what we can say is the following: 

Theorem 1.1. Truncated Pascal matrices are invertible iff the following equiva- 
lent conditions hold: 

• r <x (i.e., Vi < Xi for all i). 

• There is no zero diagonal entry. 

The primary goal of this paper is to verify this result. As we do, we develop 
some connections between univariate polynomial interpolation. Following this, we 
derive expressions for the determinant of these truncated matrices, which involves a 
connection to multivariate polynomials. The paper is organized as follows: 

1. In section (2), we use elementary matrix methods to show that r < a; is a 
necessary condition for invertibility. 

2. In section (3), we show that the invertibility of truncated pascal matrices is 
equivalent to a particular two-point polynomial interpolation problem. 

3. In section (4), we investigate the two-point polynomial interpolation problem 
in more detail. For this, we connect the problem of invertibility to the theories 
of BirkhofF interpolation and Polya systems. In doing so, we develop a new 
construction of incidence matrices using a Boolean "sum-dot" operation for 
building up and trimming well-poised incidence matrices. 

4. In section (5), we establish that invertibility conditions derived considered 
above is both necessary and sufficient. To do so, we use the machinery pre- 
viously in section (4). 

The connection between truncated Pascal matrices and polynomial interpolation 
described in sections (3) and (4) is known as BirkhofF interpolation. These latter ma- 
trices are invertible iff they satisfy a certain condition, known as the Polya condition. 
This reference goes back to a paper of G. Polya in [5], as used by J.M. Whittaker 
in [4], and later generalized by D. Ferguson in [2]. As shown, it turns out that the 
Pascal matrix is fundamentally connected to Hermite polynomial interpolation, while 
truncated Pascal matrices are fundamentally connected to 2-point Birkhoff interpo- 
lation problem studied by Polya. The more general BirkhofF interpolation problem 
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(with more than two interpolation points) was presented by G. D. Birkhoff in [T]. 

2. Necessary Conditions for Invertibility of Truncated Pascal Matri- 
ces. In this section we derive very simple necessary conditions for the invertibility 
of truncated pascal's matrix. The proof is rather straight-forward as well. Later we 
show that these conditions are sufficient as well. 

Lemma 2.1. Let A be a n x n square matrix. If A[k : n, 1 : /c] is a zero block for 
some k ^ 1 : n, then A is not invertible. 

Proof. The proof is by induction. If A is a 1 x 1 matrix with a zero on a diagonal 
element, i.e., A= [0], then it is clearly singular. The same is obviously true for 2 x 2 
matrices - here the possibilities are for fc = 1 and = 2, respectively: 



"0 *■ 








or 




* 








where * indicates any number (zero or nonzero). Now, suppose that the the result is 
true for all (n— 1) x (n— 1) matrices, and consider an n x n matrix A with A[k : n, 1 : /c] 
is a zero block. Let Mk be the minor deleting the first row and fc-th column. Then, 
Mfc is a matrix of size (n — 1) x (rt — 1) with a zero block to the left and below a 
diagonal element. Hence, expanding the determinant along the first row gives a zero 
determinant. Hence, the result is true for all n. □ 

For example, consider: 

* H< H< * 
^ H< H< H< H< 
~ * ' 
.0 *. 

Since the diagonal element A(3, 3) is zero, as well as the block to the left and below 
this element, then the following minors (expanding along the first row) each have a 
zero block below and to left of a diagonal element: 



* 


* * 










* 


* * 













* 







* 







* 


and 











.0 







.0 


*. 




.0 







.0 





0. 



Hence, the determinant of these minors (and original matrix) is zero, by induction. 

Now, we can derive a simple necessary condition for invertibility of truncated 
Pascal matrices. 

Theorem 2.2. Let r — [tq, . . . ,rd] and x = [xq, . . . ,Xd\ be selections (strictly 
increasing sequences with ro > and Xq > 0) of the rows and columns, respectively, 
of the Pascal matrix T . IfT{r, x) is invertible, then r < x. I.e., rfc < Xk for k — : d. 
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Proof. Recall that the truncated Pascal matrix can be written 

T{r,x) = 



with i ranging over the rows and j over the columns, for i, j ~ : d. Assume that 
> Xk for some k. Since the sequence r and x are strictly increasing, we have: 

r<j > Td-i > ■ ■ ■ > rk > Xk > Xk~i > ■■■ > Xq. 

Therefore, for j < fc < i we have > Xj. It follows that (^^) — for these i and j, 
which gives a lower zero block below and to the left the k-th diagonal element. By 
lemma I^TTl T(r,x) is not invertible. □ 

Noting that a diagonal element (^^) =0 is zero for some fc iff > a;^, we have 
the following corollary. 

Corollary 2.3. IfT(r,x) is invertible, then all diagonals elements are nonzero. 

In the next sections, we will prove that these conditions are sufficient for invert- 
ibility. 

3. Truncated Pascal Matrices and Polynomial Interpolation. Let A be 

the sequence of functionals (data map) 

.6^ S^D 5r,D' , 

A, 



0!' 1! ' ' z! ' ' 
with action 

and let V be the sequence of monomials 

V =[l,x,x\...]. 

We will be interested in the maps Aq and Ai, when a = and a ~ 1, respectively. 
In particular, it is straight-forward to see that T = A'^^V. It follows that, for any 
selection of rows and columns of T, 

T{r,x) = {AjV){r,x)=AiirfV{x). 

And so, invertibility of truncated Pascal matrices is equivalent to correctness in poly- 
nomial interpolation. We state this formally in the following lemma. 

Lemma 3.1. Let r = [ro, . . . , r^] and x = [xq, . . . , Xd] be selections of the rows 
and columns of T , respectively, with < < r^+i and Q < Xi < Xi^i. Then, 
T{r,x) = Ai{r)'^V{x). 
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Let X be the complement of a; in [0 : n], with n > max{r(d), x{d)}. TypicaUy, we 
choose equaUty. Now, consider the data map 

A:= [Ai(r),Ao(s)] 

and basis 

V:=v{x,x) = 1^=° , . . . , r° , . . . , t^"--^] . 

Then, we have the following: 

Lemma 3.2. T{r,x) is invertible iff A^V is invertible. 
Proof. As illustrated below, A^V is the block matrix 

rV = [Ai{r),Ao{x)fV{r,x) = 



Aj{r)V{x) 


Kl{r)V{x) 




T{r,x) 


B 


Ki,{x)V{x) 


K'^{x)V{x) 







I 



for some matrix B. Therefore, det(A-^y) = det(T(r, x)), and so A^V is invertible iff 
T{r\ x) is invertible. □ 

For example, let r = [0,2,4,7] and x = [1,2,5,8] in [0 : 8]. Then, with x = 
[0,3,4,6,7], we get 



A=[Ai(r),Ao(c)] 



5i SiD^ SiD'^ SiD^ §0 SqD^ SqD'^ SqD'^ SqD'' 
0! 2! 4! 7! 0! 3! 4! 6! 7! 



and 



V=V{x,x)=[t e 1 f6 i7j_ 

And so, KT'V becomes: 



t 




i5 








t4 


i6 














il) 


il) 


it) 


il) 


il)] 


il) 


Q) 


il) 


il) 


i^) 


il) 


it) 


il) 


il) 




il) 


il) 


il) 





il) 


it) 





il) 


Cr) 


(?) 


(?) 


(?) 


(?) 


(?) 


(^) 


(?) 


Cr) 














1 





























1 





























1 





























1 





. 























1 



The connection established in the previous lemma, will be used to determine 
necessary and sufficient conditions for the invertibility of truncated Pascal triangles by 
considering a two-point interpolation problem. We describe the interpolation problem 
in the next section. 
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4. Birkhoff Interpolation and the Polya condition. As it turns out, the 
kind of interpolation needed here to solve this problem is what is called (2-point) 
Birkhoff interpolation. In general, Birkhoff interpolation concerns when functionals of 
the form Sx^D"'^ are linearly independent on polynomial spaces. For example, we can 
interpolate a quadratic polynomial with the data /(O), /'(O) and /(I), but not with 
/(O), /'(I) and /(2). Here, we are interested in the two point Birkoff interpolation 
problem, with a:o = and xi — 1. To solve this problem (in the more general setting) 
Ferguson (PJ) used incidence matrices. 

Definition 4.1. An incidence matrix E for 2-point interpolation problems on 
n„ is a 2 X (n + 1) matrix 



E 



eoo eoi 
eio en 



eir. 



of ones and zeros, with exactly n 
j ^0 -.n, and X^Lo E"=o ^ij = " 



1 ones. I.e., Cij = I or for all i — : I and 
1. 



For example, the following are incidence matrices, each of dimension 2x6 with 
exactly 6 ones: 



1 
1 



1 
1 1 



10 11 
110 10 



Here, n = 5. As in [2^, we define Mj to be the cumulative column sum 



i=0 fc=0 



for j = : n. In the previous examples, M = [1, 2, 3, 5, 5, 6] and M = [0, 2, 3, 3, 5, 6]. 
Note that M5 = 6 is the total number of ones in both examples. In general, M„ = n+1 
for any incidence matrix. The problem of 2-point Birkhoff interpolation problems was 
studied by Polya. With respect to the incidence matrices (of any size), the following 
defintion is used. 

Definition 4.2 ([2])- The incidence matrix E satisfies the Polya condition if 
Mj > j for j ^ : m. 

So, in the above two examples, the first matrix is Polya, however the second is not 
because Mq = < 1 and M3 = 3 < 4. The following result, proved independently by 
Polya and Whittaker (as also described in [5]), gives necessary and sufficient conditions 
for correct interpolation. 

Theorem 4.3 (adapted from [3], [3]). Let E be a 2 x {n + 1) incidence matrix 
with entries etj, ones or zeros. Let 

A [SuD^ : e,, = 1] 
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and V" = . . . Then, the system A^V" is invertible iff E satisfies the Pdlya 
condition. 

We now establish some results concerning the boolean algebra of incidence ma- 
trices. The Boolean operations of multiplication (intersection), addition (union) and 
conjugate (complement) are defined by 

1-1 = 1 and 0-0 = 0- l = l- = 0, 
+ = and 1 + = + 1 = 1 + 1 = 1, 
= 1 and T = 0, 

respectively. Then, for Boolean matrices Ei and E2, we define 

E\\E2 



E^{1,:) + E2{1,:) 
.^1(2,:) -^2(2,:). 



For example. 



10 10 
10 110 1 



10 11 
110 10 



10 111' 
1 



The identity / such that E = E+I = I+E for all 2 x (n + 1) Boolean matrices is: 





1 1 



Moreover, for any Boolean matrix E — E+E. 





1 1 



Lemma 4.4. Let Ei and E2 be two 2 x (n + 1) incidence matrices. Then E 
E1+E2 is an incidence matrix iff 



\E^{l,:)-E2{l,:)\ = \Ei{2,:)-E2{2,:)\. 
Moreover, E is Pdlya if both Ei and E2 are Pdlya. 

Proof. Since Ei is an incidence matrices, we have (n + 1) = \Ei\ 
Ei{l, :)| + \Ei{2, :)|, and likewise for E2. Then, 

\E\ = 1^1(1, 01 + \E2{1, 01 - |£;i(l, • E2{1, 01 + (n + 1) 



(|i;i(2,0| + |i;2(2,0|-|i^i(2,0-^2(2,OI) 



)l + 1^^2(1, 

)l + |£^2(l, 

)\ + \Ei{2, 



)| + (n + 1) - {\E[^\ + I^^IMI) iff hypothesis 

)| + (n + 1) - ((n + 1) - |i?i(2, 01 + (n + 1) - \E2{2, :)\) 
)\ + \E2{l,:)\ + \E2i2,:)\-{n+l) 



(n + 1) + (n + 1) - (n + 1) = n + 1. 
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For the second part, we first note that if Ei and E2 arc both Polya, then with 



Mo = Ml + M2 = \elo + elo\ + |e}o • e?o| 

either one of ego is 1, or otherwise both of e^g are 1 since Ei and E2 are Polya. 
Proceeding inductively... □ 

The converse to the last statement is not true. That is, E Polya does not neces- 
sarily imply that Ei and E2 are Polya. For example, consider 



E 



110 10 

10 10 

1 1 0' 
1110 



110 

10 10 1 

1 0' 
10 111 



10 1 
10 11 



In the first decomposition of E, the two matrices are Polya, and in the second they are 
not. Therefore, there are at least some Polya matrices that can be broken down into 
a sum-dot of two incidence matrices, and these may or may not be Polya. But the 
question here can one always express incidence matrices as the sum of two incidence 
matrices with less ones in the first row, and can one do it with Polya matrices. We 
will answer this question with a constructive proof. Suppose that E is an incidence 
matrix, and suppose also that |£^(1, 01^2 (i.e., there are at least two ones in the first 
row. Then define Ei and E2 as follows: 



Construction: Let d 
except the following: 



|£(1,:)|. Take Ei and E2 equal to E for all elements 



1. For El, exchange the first one in row one of E with the first zero in row two. 

2. For E2, exchange the last d—1 ones in row one of E with the last d—1 zeros 
in row two. 



In the extreme case, when d 
examples: 



0, we get El = E2 = I- Consider two other 












1 








0' 

























+ 











1 








0" 


1 


1 


1 





1 


1 


1 




1 


1 


1 


T 


1 


1 


1 




1 


1 


1 





1 


1 


1 



and 



1 1 1 0' 
1110 10 



OJ 1 1 

1 1 1 |T| 1 



10 
1110 



o' 

1 



All matrices are Polya. The entries that were exchanged from the original E are 
indicated in boxes. In the first example, since d = 1 we get Ei = I and E2 = E (i.e.. 



there are d—l = entries exchanged in the second). Note also that :)| = \Ei{l,: 
)| + \E2{1, :)|, which will be useful to us later. In general, we can say the following: 

Proposition 4.5. Let E he an incidence matrix. Suppose that E\ and E2 are 
constructed as above. Then, 

1. E = EitE2 

2. El and E2 are incidence matrices. 

3. If E is Pdlya, then Ex and E2 are Pdlya. 
I \ E{1,:) \ = \ Ei{l,:) \ + 

5. \E{2,:)\ = \Ei{2,:)\ + \E2{2,:)\ 

6. Ifm = 0, E = El = E2 = I. 

7. Ifm=l,Ei=I and E2 = E. 

8. Ifm>l,E,Ei and E2 are distinct. 

Theorem 4.6. Every incidence (Pdlya) matrix can be written as a sum-dot of two 
incidence (Pdlya) matrices. If m>2, the sum can be non-trivial (E, Ei and E2 are 
different matrices), and one can have |£Ji(l,:)| < |£J(1,:)| and |£'2(1,:)| < |£^(1,:)|. 

Let X be the complement of a; in [0 : n]. That is, if x = [2, 4, 5] with n — Q, then 
X = [0, 1, 3, 6]. We define E^.x the 2 x n + 1 matrix with £^^^(0, j) = 1 if = 1, and 
i?r,a(2, j) = 1 if = 1, with all other entries 0. 

Lemma 4.7. Let Er^x be a 2 x (n + 1) matrix as defined above. Then, it is an 
incidence matrix. 

Proof. We need to show that there are exactly n + 1 ones (hence, n + 1 zeros) . 
We denote this as |-Br,cc|. But this is easy, since 

\Er,x\ = \r\ + \x\ = \r\ + {n + l)-\c\ = {d + l) + {n + l)-{d+l)=n + l. 

□ 

Next, we want to consider combining incidence matrices in certain ways, i.e., by 
addition. That is 

5. Sufficient Conditions for Invertibility of Truncated Pascal Matrices. 

We are now ready to determine conditions for the invertibility of truncated Pascal 
matrices. We provide two proofs, both involving Polya systems. The first using an 
inductive argument along with the sum-dot decomposition described in the previous 
section, and the second is longer but more direct, using a counting argument. 

Theorem 5.1. Let r = [rg, . . . , r^] and x = [xq, . . . , x^] he selections (increasing, 
1-1 maps [0 : d] ^ [0 : n]) of the rows and columns, respectively, of the Pascal matrix 
T, for some d < n. The matrix A := T{r,x) is invertible iff r < x (i.e., ri < Xi for 
i = 0:d). 

9 



Proof. Proof 1: Fix n. Assume by way of induction that the resuh is true up 
to c? — 1 (i.e., + d — 1). Let tq, . . . ,rd and xq, . . . ,Xd as described in the theorem. 
Define the matrix with first row E{1, r{i)) = 1 for i = : d, aU other entries zero, and 
second row E{2, x{i)) — 0, all other entries 1. Then, there are d+l + {n — d) = n + 1 
ones. Hence, this is an incidence matrix. The first row corresponds to 5iD^ , and the 
second row 5qD^ . Now, define incidence matrices Ei and E2 as in the construction 
given earlier, so that E = E1+E2, with — d and |i?2(l,:)l ^ 1- Assume 

that Tfc < Xk for all k. Then, the same is true of the submatrices. By the inductive 
hypothesis, Ei and E2 are Polya. By Lemma 14.41 E is Polya. Therefore A„14, is 
invertible iff r < a;. By lemma |3^ T{r,x) is also invertible iff r < a; 

Proof 2: Our strategy in proving this theorem is to expand T(r, x) to an {n + 
1) X (n + 1) matrix B by adding diagonal entries where we want the coefficient of P 
to be zero. Then, we show that B{r,x) = T(r,x), and that B is invertible iff B{r,x) 
is invertible. 

Let X be the complement of x in [0 : n] , a vector of length n — d, and let E be 
the 2 X (rt + 1) incidence matrix with E{0,Xj) = 1 ioi j — : n — d and £'(1, rj) — 1 
for j = : d, with all other entries 0. For example, suppose that n — 6, x = [1, 3, 5] 
and r = [1, 5, 6]. Then with x ~ [0, 2, 4, 6], we get 



E 



10 10 10 1 
1 1 1 



In particular, since #r = #c, the total number of ones is 



M„ = #r + #x 



(d+l) + (n+l)-(d+l) = n+l 



as required. Now, we have 




j=0 



and 
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Therefore, if rk > Xk for some k, then we have 



Ck rk- 1 

= {xk -k) + {rk-i -Xk) + k 
= rk-i- 

However, the Polya condition requires that M^^-i > rk- Hence, the system is not 
Polya when rk > Xk for some k. 

Now, suppose that < Xk for all k. Let j G [0 : n]. Then, let ko be the highest 
k such that Xk < j- Then, r^o < Xko < j- And so, 

j j 

Ik 



fc=0 fc=0 



fc=0 fe=0 

Xka ^(=0 

= Y^eok + (j - a^fco ) + XI ^1*^ 
fe=o fc=o 

= (a;fco - ko) + {j -Xko) + {ko + 1) 

= j + l. 

Prom this we conclude that E is Polya iff rk < Xk for all k. 
Now, let 

A„ := [(5i£>-'' : Cij = 1]. 
By the above theorem, the system 

Hn = AlV" 

is invcrtiblc iff E satisfies the Polya condition. That is, wc can say = =^ a = 
iff is Polya. But if A, = SqD^ for some functional in A, this implies that aj = 0. 
Such coefficients correspond to the first row of E where there ones, i.e., those entries 
in X, Therefore, we conclude that Hna = a = iff Hn[r,x]l3 = => /3 = 0. 
Therefore, Hn is invertible iff Hn[r,x] is invertible. Since Hn = TnDn, we get that 
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T(r, x) is invertible iff Hn is invertible, which is true iff E is Polya, which is true iff 
rfe < Xk for all k. This proves the claim. □ 

Corollary 5.2. The matrix T{r,x) is invertible iff all diagonal elements are 
nonzero. 
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